Framework for discrete-time quantum walks and a symmetric walk on a binary tree 
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We formulate a framework for discrete-time quantum walks, motivated by classical random walks 
with memory. We present a specific representation of the classical walk with memory 2 on which 
this is based. The framework has no need for coin spaces, it imposes no constraints on the evolution 
operator other than unitarity, and is unifying of other approaches. As an example we construct a 
symmetric discrete-time quantum walk on the semi-infinite binary tree. The generating function 
of the amplitude at the root is computed in closed-form, as a function of time and the initial level 
n in the tree, and we find the asymptotic and a full numerical solution for the amplitude. It 
exhibits a sharp interference peak and a power law tail, as opposed to the exponentially decaying 
tail of a broadly peaked distribution of the classical symmetric random walk on a binary tree. The 
probability peak is orders of magnitude larger than it is for the classical walk (already at small n). 
The quantum walk shows a polynomial algorithmic speedup in n over the classical walk, which we 
conjecture to be of the order 2/3, based on strong trends in data. 

PACS numbers: 03.67.Ac, 89.70.Eg, 02.50.Ga, 05.40.Fb 



I. INTRODUCTION 

Random walks on graphs (Markov chains) are used 
extensively in science. They provide a number of now 
standard approaches and models in physics. Application 
of such ideas to evolution of quantum systems has led 
to the emergence of the field of quantum walks, princi- 
pally distinguished between discrete-time (DTQW) and 
continuous-time (CTQW) quantum walks. However, be- 
ing unitary (reversible) processes, quantum walks are 
very difi^erent from their classical stochastic (Markovian) 
counterparts. 

Quantum walks are used to approach varied problems, 
for instance, on quantum lattice gases, arrow of time, 
generalized quantum theory, exciton trapping, or topo- 
logical phases [1-5] . They may become a general tool for 
building physical models; for example, see a summary in 
[5] and CTQW in transport phenomena [6]. In quantum 
computing, they are a universal primitive [7], and in the 
algorithmic context, the principle alternative to quantum 
Fourier transform. The field has developed since its initi- 
ation [8-12], with established algorithmic uses, examples 
of dramatic superiority over classical approaches [13, 14], 
and implementations [15]. See [16, 17] for a review and 
a recent summary. 

Standard approaches to quantum walks generally stem 
from memoryless classical walks. However, since quan- 
tum evolution is memoried (unitary), it seems natural to 
approach construction of quantum walks from classical 
random walks with memory. Relation between unitarity 
and memory in walks has been noted [1, 3]. (Also, in 
computer science memoried and biased approaches are 
common and beneficial algorithmically.) 
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It is our observation that DTQW are most directly 
related to classical walks with memory 2. In this paper 
we present a general DTQW framework as a direct analog 
of a specific representation of memory-2 classical walks. 
With it we construct a symmetric DTQW on a binary 
tree, starting from a pure state at an arbitrary level in 
the tree, and compute its amplitude at the root. 

Sec. II starts with a representation for memoried clas- 
sical walks, that is particularly interesting for quantum 
walks due to the specific form of the Markov tensor. In an 
analogy with it, we then define a framework for DTQW. 
Walks are built by choosing the evolution operator with 
no constraints other than unitarity. They evolve in the 
product of state spaces, while the key component of the 
operator acts on single states. There is no need for "coin" 
degrees of freedom. The framework is flexible and suit- 
able for general graphs. It is also unifying of other ap- 
proaches, notably coined and Szegedy's [18]. 

In Sec. HI we apply this framework to a binary tree, 
a structure with many uses in physics. It is a natural 
environment for quantum walks, but difficult to utilize 
with current techniques for DTQW. A successful specific 
construction exists for CTQW [13]. 

We construct a symmetric walk on the semi-infinite bi- 
nary tree, and calculate its amplitude at the root, as a 
function of time and of the initial level in the tree. This 
involves path enumeration using regeneration structures, 
manipulated with the z-transform. The obtained closed- 
form generating function yields the analytic asymptotic 
for the amplitude, which we also compute numerically. 
The amplitude has a sharp peak and a power law tail, 
completely unlike the corresponding classical random 
walk, and shows a polynomial speedup in n. The data 
strongly suggests the order of the speedup of 2/3. 

A summary is in Sec. IV, and appendices discuss the 
steepest descent calculation, and the classical walk. 



II. A FRAMEWORK FOR DISCRETE- TIME 
QUANTUM WALKS 

The main approach to DTQW fohows ideas of classical 
memory less walks, and needs an auxiliary "coin" degree 
of freedom. Here we employ a specific representation 
of memoried Markov chains (11 A), as an analogy for a 
general DTQW framework (II B) which does not need 
coin spaces. Relation to other approaches is discussed at 
the close of this section. We start with an example that 
serves as a motivation for using memoried walks. 

Memoried walks and coined DTQW. Classical walks 
with memory are walks with internal states: apart from 
its state, the walk carries other information. Here we 
use one such walk in one dimension, with the following 
property: the next step depends on the direction of the 
previous one. If the walker came to a site i from the 
site i — 1, the probability to go to i + 1 (to maintain the 
direction) is p, while the probability to go to i— 1 (reverse 
direction) is 1 — p. This is often called a persistent walk. 
We now show that coined DTQW are a special case of a 
quantum analog of classical persistent walks. 

Consider a standard coined walk on a d-regular graph 
with n vertices (for example, [12]). The state of the walk 
is in the direct product of two Hilbert spaces: an auxil- 
iary ("coin") space Ha, spanned by d states \a), in which 
a unitary operator C mixes components; and a space 
of vertices, Hv, spanned by n states \v). The evolu- 
tion operator acts in this product space Ha ® Hy as: 
U (|a) (g) \v)) = S" (C (g) /) (|a) ig) \v)). On a cycle with the 
Ha basis {|t)Ji)}7 the shift S can be implemented as 

s = It) (tl ® Ej \j + 1) {j\ + \i) ai ® E, |j - 1) (il- Now 

consider for C a generalized Hadamard coin operator, 



C = 



\/i - p -Vp . 



and evolve the state. Starting from the pure "up" state, 

S-iC®I)\t)®\i) 

= ^ • (Vp It) ® K) + v/T^ It) ® N) 



= Vp It) ® K + 1) + yr^ It) ® N - 1). 

For the pure "down" initial state. 



(1) 



S-{C(g>I)\i)(g,\i) 

= v/r^it)®i* + i)-Vpit)®N-i). (2) 

This is a persistent walk: it maintains the direction with 
probability p, and changes it with probability 1 —p. The 
obtained walk undergoes the same spectral analysis as 
its unbiased special case {p = 1/2, the Hadamard walk). 
Walks with a general coin have been studied (for exam- 
ple, [19-22]), as well as persistent (correlated) walks and 
their relation to DTQW [23-25] . With this example we 
point out the direct correspondence between them. Note 
that the directionality of the walk shows up as soon as 



the coin transformation is allowed to have p ^ 0.5. In 
other words, the standard Hadamard transform generally 
implements a persistent walk (rather than a memoryless 
one), only with equal probabilities. 

A. A representation for classical memory— 2 w^alks 

Walks with memory 2 are such Markov processes where 
the next step depends on two states: the current one, 
and the previous one. Walks with memory are generally 
studied by using a suitably enlarged state space. In par- 
ticular, a memory-2 Markov chain can be represented as 
a memoryless one over the space with n^ states. The 
transition matrix is then large (n^ x n^) and sparse. 

Instead, here we represent a Markov chain with mem- 
ory fc by a probability distribution ^{t) of dimension k, 
while the Markov tensor M is then of dimension k + 1. 
For a memory-2 walk over n sites, the space has dimen- 
sion n, and each state is labeled by two indices (the site 
the walker came from, and the current site). So the prob- 
ability distribution is two-dimensional. 



fi{t) 



Mo,o(0 



l^0,n-l{t) 



Mn-l,o(i) ■■■ Mn-l,n-l(^) 



The matrix fj,ij can also be given by a column of rows 
ri, or by a row of columns Cj, which we use below. The 
following representation for the third-rank tensor Ai and 
its action is convenient. 

Let Pij\k be the conditional probability for the transi- 
tion j —>■ k, given that the walk came to j from i. All 
transition probabilities {pij\k} define the evolution oper- 
ator A4 = [PqPi ■ ■ ■ Pn-i], as n layers of n x n transition 



matrices Pj, j — 0,1, 



1, one for each site: 



P,= 



POjIO P0,i|l 



Pn-l,j\0 Pn-l,j\l 



P0,j\n-1 



Pn — l.j\n — l 



Pj are by construction transition probability matrices, 
and this is the only requirement imposed on them. 

The evolution of the state, /it+i = /it A^, with M act- 
ing to the left, is defined as 

^l{t)^^i{t + l) : rj{t + l)^c]{t)P„ 

for each j = 0, 1, . . . , n — 1, where rj and cj are the j-th 
row and transposed column, respectively, of the matrix 
/i. In words: At each site j , the Pj associated with that 
site acts on the transposed j-th column of ^i{t), giving 
the j~th row of the evolved fi(t + 1). 

Instead of an n^ x n^ probability matrix, we use n of 
nxn probability matrices Pj . They implement the evolu- 
tion: the j'-th column oi fj,(t) has probabilities to arrive to 
j from any site, and after the action of Pj the j-th row of 
fi{t + 1) has probabilities to go from j to any site. Thus 



action of all transition matrices on all columns evolves 
the probability distribution over all paths. The stochas- 
tic nature of the process is carried by the assignment of 
{Pij\k} transition probabilities in Pj matrices [26]. 

The Pj transition matrices are simple in most cases of 
interest. Consider the cycle, a space {0,1,..., n} with 
identified ends (0 and n), with only nearest-neighbor 
transitions, {j ± l,j) — > (j,j ± 1). Take the persistent 
walk, with probability p to continue, and 1 — p to reverse, 

(i - 1, j) ^ {j,j + 1), with probability p 

{j - 1, j) -^ {j,j - 1), with probability 1 - p. 

To obtain this walk, the Pj matrices have the following 
block centered at {j,j) ( mod n), 



Pj = 



1-p 




P 



P 


1-p 



(3) 



The rest of the diagonal of Pj has Is, other elements are 
0, except for the transitions between sites = n, and n— 1 
or 1 (boundary conditions), which are Po,n-i|o — ^ ^ P 
(reverse), Po,n-i|i = P (continue), etc. Action of these 
Pj by the above prescription carries the walk. 



B. Quantum walks: the interchange framework 

The above classical procedure for memory-2 walks is 
directly elevated to define quantum processes. 

Consider a basis in an A^-dimensional Hilbert space, 
with vectors labeled as {\i),i = 1,2,...A^}. They rep- 
resent states that the walk is performed on, enumerated 
on a general graph. The state of the walk is given in the 
product C^ X C^ spanned by these bases, by states at 
the previous (|i)) and current (|j)) step, 



m)) = Y,c^^j{t)\i)®\j). 



(4) 



The evolution is specified by 

m + l))=UXm)), \yj{t)}^{UXY\^{0)}, (5) 

where X is the interchange operator and U is defined via 
unitary operators Uj in C^, assigned for each site, 

X: \i) ® \j) ^ \j) ^ \i) 

N~l 



j^ivj, where H^ = |j)(j|. (6) 



IIj selects the first state in the product, and Uj acts on 
the second. Before explicit examples, we make a few 
general comments. 

Consider a pure state of the walk \i) (g) |j), represented 
by an arrow pointing from the previous («) to the current 
(j) site. The interchange initiates the walk forward by 
"reversing the arrow." Then the Uj operator distributes 
the "tip of the arrow" to all sites the process can ac- 
cess (generally in the subspace of adjacent nodes), and 
the evolved superposition is obtained. This is best seen 
in the forthcoming example of the binary tree (Fig. 2). 
The explicit reversal X is crucial; then Uj completely 
controls the evolution over site j, by acting on the origi- 
nating state (s) |i), and sending the process over all paths 
to a new state. The framework does not place any condi- 
tions on these operators, except for unitarity of quantum 
evolution. We are free to choose (or construct) them as 
needed to implement quantum walks. 

Note that this construction needs no mention of clas- 
sical processes. The representation of classical memoried 
walks in Sec. II A is given for motivation and insight, 
and we now comment on this relation. The discussed 
interplay between interchange and (local) Uj, critical 
for this formulation, has a clear analog in the classical 
representation — recall the transposition before (local) Pj 
evolve the distribution. Also, the freedom to craft any 
(unitary) Uj to implement quantum walks corresponds 
to a classical property, as Pj may be any (probabilistic) 
matrices. Finally, the classical representation has no ex- 
plicit coin toss, and there is no need in the quantum case 
to mimic randomization via a coin degree of freedom; 
here Uj drive the walk and mix components. 

Relation to other approaches. For a comparison 
with a memoryless (coined) approach, consider a walk 
on the line, over the state space S = {\j),j = 0, 1, . . .}. 
The quantum walk (1) and (2) is obtained with 



Ujip) 
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Vp 
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Vi-p 


1 



J=0 



with the block centered at (j, j). The rest of the diago- 
nal has Is, and other elements are 0. (On a cycle there 
are elements needed for boundary conditions, like in the 
classical case.) The square roots provide for probability 
being the square of the amplitude, and the —y/p sign is 
needed for unitarity. In this simple case, the choice of 
Uj follows from the classical memoried walk (3), but in 
general a classical analog is not needed. 

Now look at the evolution steps from pure states. 
Starting from \i — 1) (E) \i) (the system is in the state 



\i), having been in \i — 1) at the previous step), 

UX\i-l)(8)\i) = U\i)(g)\i-l) 

= {j2\^){j\®Uj)\z}^\z-l} 



= ^/T^\i)®\i-l)+y/p\t)(E>\i + l). 
Similarly, for the initial \i + I) (E)\i) state. 



UX\i + l)(g>\i) ^-^\i)(g>\i~^l) + ,/l-p\i)(g)\i + l). 

This is an isomorphism of the memoryless-based walk of 
Eqs. (1) and (2), via identification 

\i - 1) (S) \i) <^ It) <S) \i) and \i + I) (S) \i) ^ \i) <» \i) . 

This analysis applies to arbitrary mixed states, as each 
component is evolved separately. The choice of p = 1/2 
restores the Hadamard walk. Thus the interchange 
framework reproduces coined quantum walks on the line 
directly, with the above choice for Uj. 

Memory in quantum walks is mentioned in literature. 
For example, it was noted in study of the classical limit 
via decoherence and multiple coins [27] , and a direct rela- 
tion between coined walks and classical memoried walks 
was observed [28] . Recently a particular "quantum walk 
with memory" [29] was studied. 

An important approach directly resorting to ideas of 
memoried walks is the Szegedy walk [18], which is the 
most prominent tool in DTQW not using coin degree 
of freedom [17]. Its construction starts from a classi- 
cal Markov chain, and the resulting evolution operator 
explicitly carries classical transition probabilities. It is 
contained in the interchange framework via the specific 
choice 



(U, 



'2\/p{J:ik)p{j,ini) - 6k 



where p{i,j) need to be classical transition probabilities. 
The present approach does not require a specific form 
of the evolution operator. It is fully defined by (4)-(6) 
alone, without reference to classical walks, and quantum 
processes with desired properties are set by choosing Uj 
without constraints. Some benefits of this are seen in the 
next section, where we construct a symmetric DTQW on 
a binary tree. A Szegedy walk on a binary tree cannot 
be obtained with equal probabilities for each branch, as 
there is no (real) solution for probabilities p{i, j) such as 
to yield a 1/3 probability for the quantum walk. 

The Szegedy walk is also a translation of a memoryless 
walk into a walk with memory. The interchange frame- 
work is a direct analog of an explicit representation for 
mcniory-2 Markov chains. This is reflected in some of its 
properties, discussed above. 

The interferometry-motivated [30] scattering walk [31] 
is performed on graph edges, scattering off of vertices (or 
subgraphs). The more general formulation [32] can be 
formally reconciled with the present framework, while 



their designs and interpretations are different, and com- 
plementary. The scattering walk has recently been used, 
in its form reflecting the physical (scattering) approach, 
for certain search problems [32-34] , hinting at benefits of 
coinless algorithms. This approach has not been related 
to stochastic processes, which underlie the motivation, 
construction, and expected uses of the present frame- 
work. We showed here how a very general formulation of 
DTQW, unifying of other similar approaches [35] , arises 
naturally from Markov chains with memory. 

Before moving to a full example we remark on the ver- 
satility of this approach. It reproduces coined walks on 
the line and Szegedy's walk, and handles a walk on a 
binary tree (next section) — each by a simple choice for 
Uj. This demonstrates flexibility, and it seems that the 
framework can help with problems that so far have been 
prohibitively difficult. 



III. A SYMMETRIC DISCRETE-TIME 
QUANTUM WALK ON A BINARY TREE 

The binary tree is a common model in physics, and 
a structure of interest in quantum computing [36] . One 
of the initiating works [8] used it as a model for deci- 
sion trees, and one of the most successful algorithms [13] 
solves a particular problem on connected binary trees, 
both using CTQW. We have not seen such progress in 
using DTQW on the binary tree, even though this would 
be beneficial for many problems. This seems to be due 
mostly to trouble in handling coin spaces that are nec- 
essary for (coined) DTQW. In this section we use the 
established framework to set and calculate a symmetric 
DTQW on the semi-infinite binary tree. We orient our 
tree with the root (single starting node) at the left with 
the tree spanning to the right (Fig. 1). 






FIG. 1. Conventions used for our binary tree. 

We focus on the following basic question. The walk is 
started from a pure state at a site in the tree at a level n, 
and we examine its amplitude at the root as a function 
of time (step) and the initial distance n. 

The desired symmetric walk has equal probability to 
step to either of the connecting nodes, having come from 
either direction. (The analysis remains unchanged for 
different choices of local Uj.) The state of the walk is 
given in the direct product of spaces, each spanned by 
states defined at nodes S = {\i)}. Label a node in the 
tree as j, and the nodes connected to it as ii (to its left, 
toward the root), i2, and 13 (to its right, away from the 




root), as in Fig. 2. Consider an evolution step from a 
pure state at j, for example, 1-00) = 1*2) '9 \j)- The action 
of the interchange X reverses the state. Next we want to 
write down the Uj matrix, acting on |i2), such that the 
evolved state has equal probabilities for either branch. 
Formally Uj operate in the space of all nodes, but they 
are reduced to the subspace of nodes to which transi- 
tions are allowed; here the adjacent ones. Thus Uj can 
be written in a block-diagonal form, with the non-trivial 
transition matrix U'^°'^ in {|ii), |i2), lis)}, and an iden- 
tity matrix over the remaining dimensions. The matrix 
elements need to satisfy the unitarity of Uj and equal 
squared amplitudes of components of the state evolved 
by it. The obtained evolution operator, with (reduced) 
transition matrix U'''°'^ in the basis {|ii), 1*2), Na)}; is 



U, 



ur 







where a 



I 



o27ri/3 



with U'""^ = ^ 

' \/3 



"l 


a a 


a 


1 a 


a 


a 1 



, (7) 



This representation holds for graphs 



of any degree, where dimensions of U and I change. 
At the root the walk can only get reflected, which is per- 
formed by interchange X; then Uq is the identity ma- 
trix. This will be accounted for. For all other states, we 
now follow the prescription of Eqs. (5) and (6). With 
U — J2ieS ^i^Ui, the step is 



I Vi) = UX l^o) =lj2\^){^\®uAx |Z2) 



^\j} 



li)®t/,N2>- \j) 



V3 



(a,l,a,0,...)T. (8) 



Thus the state is evolved by UX to the superposition 



1*2) (^ Ij) ^ Ij) » f ^|ii) + ^|Z2; 



V7!'"^"7! 



V3 



1*3; 



(9) 



Each component of the superposition takes the next step 
from its node in the same way, and the process spreads 
over the tree [37]. 





FIG. 2. A step taken from a pure state \i2)<^\j) [Eqs. (7)-(9)]. 
The state is sent by UX over all available paths. Probabilities 
for either branch are chosen to be equal, regardless of how 
the walk approaches the site j (the walk is symmetric). Each 
component of a general (mixed) state is evolved this way. 

In comparison with Markov processes, a quantum walk 
can be considered as the evolution of the amplitude dis- 
tribution. Also, the causality typical of the local dynam- 
ics of the classical Markov evolution, seen in Sec. II A, is 



reflected in the quantum walk. Note how the concerted 
action of X and Uj implements the "arrowed" (memo- 
ried) nature of the evolution, mentioned in Sec. II B. 

For organizing the calculation, it is useful to note 
the connection between directionality and weights of the 
components of the evolved state. The component that re- 
verses the direction of the previous step has the coefficient 
l/\/3, while the other two have a/\/3 (see Fig. 2). This 
is always the case for this walk (not only for [22) <8) |j)), 
since it is symmetric, as explicit in Eqs. (7)-(9). 

Outline of the calculation. The amplitude at the root 
at time t is computed as a sum of the contributions (am- 
plitudes) of all possible (classical) paths that are at the 
root at that time. This is practically a discrete form of 
the path integral in quantum mechanics, and is a stan- 
dard technique [1, 11]. So we count all such classical 
paths on this structure, weighted appropriately. 

The presence of a reflective boundary (the root) com- 
plicates the classiflcation and counting, and we use re- 
generation structures, which are then handled via the 
z transform. The obtained explicit expression for the 
transform is complicated, and analytically the asymp- 
totic of its inverse is found, using the method of steepest 
descent. The full amplitude is calculated numerically. 

For brevity in involved descriptions, we sometimes use 
"paths h{t)" to refer to "those paths that contribute to 
the part of the amplitude (that is named) h{t)P 



A. Path counting and regeneration sums 

Enumeration of paths, weighted with appropriate co- 
efficients (amplitude) , is a combinatorial problem. Given 
the symmetry between up and down directions, the tree 
can be projected to a line bisecting it. The paths on 
the tree can be classified, and this results in rules for an 
equivalent walk on that line. 

A component of the state at a site is directed either 
toward or away from the root; and it can either continue 
in the same direction or reverse it in the next step. For 
example, the component directed toward the root (to the 
left) can continue toward the root (taking the branch to 
the left), with the amplitude a/v3, or it can turn and 
step away from the root, by directly reversing or (and) 
by taking the other branch leading away, with the total 
amplitude of (1 + a)/y/i [see Fig. 2 or Eq. (9)]. 

Summarized by the direction of the previous step, the 
walk on the line can take the next step as follows: 
When directed away from the root (to the right), it can: 

(i) turn back, with the coefficient 1/^/i (left turn); 

(ii) continue, with (a -I- a)/^/2> (right step). 
When directed toward the root (to the left), it can: 

(i) turn away, with (1 -I- a)/-\/3 (right turn); 

(n) continue, with c/a/S (left step). 

There is a special case, not following the above clas- 
sification, which complicates the counting of paths con- 
siderably. We need to count weighted paths that are at 



the root at time t. Paths generally reach the root in 
fewer than t steps, then going back and forth in the tree, 
possibly touching the root again in the process, before fi- 
nally finding themselves at the root at time t. Whenever 
they touch the root their next step can only be a turn 
back, with the coefficient 1, and this docs not fall into 
the above classification. To account for it the paths need 
be enumerated particularly carefully. 

All paths that are at the root at time t have the follow- 
ing structure. They touch the root for the first time at 
one point (step s), and we call the amplitude for this part 
of the path hn{s). Then they go out in the tree, eventu- 
ally coming back to the root at step t, possibly touching 
it multiple times in the process; we call the amplitude of 
this part of the path G{t — s). This is encoded by the 
convolution over the first contact with the root, and the 
amplitude, represented by weighted paths that are at the 
root at step t, starting from level n, is 



counting. We work with them using the z transform, 

oo 
t=0 

Applying this transformation to Eqs. (10) and (11), 

H„{z) ^ hn{z) G{z), and G{z) = g{z) G{z) + 1, 
1 



Hn{z) = hn{z) 



1 - 9{z) ' 



(12) 



We need transforms of amplitudes of paths reaching the 
root for the first time [/i„(z)], and of simple loops \g{z)]. 
At this point we note a known combinatorial result. 
The number of paths on a lattice in two dimensions, going 
from (0, 0) to (2ri, 0), taking only Northeast or Southeast 
steps, with k peaks, is given by Narayana numbers [38], 



ff„(i) = E/ln(s)G(t-s). 



(10) 



^(^'^)^^U 



n 
fc-1 



(13) 



After the root is touched for the first time, the remain- 
der of the walk is a root-to-root path, considered in- 
dependently as G{t) [accounting for the n — Q case, 
G{t) — i?o(0]- It consists of: a "simple loop" 17(5), that 
goes from the root into the tree and back to it (reaching 
it again for the first time), followed by the rest of the 
path G{t — s), which may touch the root multiple times, 
so again comprised of simple root-to-root loops (Fig. 3). 
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FIG. 3. The number of paths is a convolution w.r.t. the first 
contact with root, with the regeneration structure of Eqs. 
(10) and (11). 



This is a convolution too, with one adjustment. To cor- 
rectly account for i = 0, it must be that 17(0) = 0, as well 
as G(0) = 1 (since t ^ ^ s = t). Thus 6o{t) need be 
added. The amplitude of such root-to-root multi-loops 

G{t) is 



G{t)^J29is)G{t's) + 5o{t). 



(11) 



s=0 



The regeneration sums, Eqs. (10) and (11), organize path 



This expression applies to the number of paths compris- 
ing the simple loops g, where peaks are positions furthest 
from the root. We first need to identify and enumerate 
"steps" and "turns" in such paths, so that we can assign 
weights to them accordingly. 

A simple loop must take an even number of steps. The 
first step is a reversal: it starts at the root, having arrived 
to it from the first node, and it can only step back onto 
the first node, so the coefficient for this step is 1. It is 
straightforward to establish that paths with k peaks take 
k left turns and fc — 1 right turns. Also, loops of t steps 
must take ^^ — {k—l) right, as well as left, steps. Loops 
with t — 2 are different: they can only step away from 
the root and return to it in the next step (left turn) ; their 
coefficient is 1 x l/VS. Thus a simple loop with k peaks, 
for i > 4 steps, bears the coefficient: 



1 \ *: /I , \ fc-1 / \ t/2-k / n \ t/2-k 

1 \ / 1 + a\ /a\ / 2a ^ 



V3/ \ VS J \V3 



Vs 



1 / l + a 



(V3)*-i V 2a2 
1 



fe-i 



(2a 



2Nt/2-l 



(V3)*-i V 2 



k—l 



(as a = e'^''^). 



For a — e^^^^^ we have 1 + a + a"^ = 0, used above. 
Summed over all possible numbers of peaks fc, and with 
the t = 2 case added, the amplitude of a simple loop is 



1 



g{t)^—5o{t-2) 



E 



(2a^Y^'-' 



t (V3)- V2-)''^(^'^ 



Using the Narayana numbers (13), with t = 2m + 2, 



m-l . 



fc=0 



m \ I rn 
k + ll\k 



Since we will need the transform of ^(i), it is helpful to 
write the above sum as an integral, using the identity 



Y.^o.pf 



k=0 



m \ m 



fc + 1 



27r 
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(l + ae*"^)"(l + /3e-*^)' 



-dx. 



Employing this, under the constraint a/3 — —1/2, 

'2a2^ ■" 



^("^^^Ti'^^^'^ + ^^^^aVs 



(14) 
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+ ae" + /3e" 



-dx. 



It is calculationally convenient to take the z-transform 
of g{t) at this point. Since loops take even number of 
steps, and g{z)t=o = .g(0) = 0, with t = 2m + 2, 



5(2)=^5(0^*=5(0)+5(2)z'+ Y. 9{t)^ 

t=0 i^4,6,... 

m— 1 

The transform of 5 is 1, and Eq. (14) becomes 
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Now we make use of J2'^=i x"/n = — ln(l — x), \x\ < 1, 
and at this point pick a = — /3 = l/v2, arriving at 



^^^^=Vl+^V3X '"^" 
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2a^ /I e* 
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Using uj = e^" and integrating by parts. 
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Here the Residue Theorem is used. The singularity at 
a; = is removable, while one of the two zeros of the 
denominator is inside the integration contour. Finally, 



?(^) 



V3 
2a2 



i + liazf 



';{az)^ + {azy 



(15) 



This closed-form expression analytically extends g{z) be- 
yond the disk \z\ < 1 on which it was defined. We now 
need to deal with hn{z). 

Paths from the n-th level that reach the root for the 
first time in s steps, with amplitude ft.n(s), first reach 
the level n — 1, generally going out into the tree in the 
meanwhile, then the level n— 2, and so forth until the root 
is hit. This is organized into paths dropping by one level 
closer to the root (with amplitude hi), convoluted with 
the rest of the walk, which itself is comprised of paths 
getting closer to the root by one level, ft,„ ~ hi * /i„_i = 
. . . ^ hi ^ . . . ^ hi (n times). Then the transform is 



hn{z) ^ hi{z) 



(16) 



Paths /ii, that for the first time reach one level closer to 
the root, are combinatorially equivalent to the paths that 
start at level 1 and touch the root for the first time. 

Consider such a path, starting at level 1 and finding 
its way to the root for the first time, in more than one 
step (Fig. 4). It must have come to level 1 from level 2, 
and it first steps back to level 2. For comparison, now 
recall a root-to-root simple loop g{s). It differs from hi 
by: the first step of g (with coefficient 1) is not taken by 
hi (which is already at level 1), and the next step of g 
(for paths with t > 2) is a right step, while the hi path 
takes a right turn. 



any (I ^ i-l) path 
hi paths 





l-l 



FIG. 4. Combinatorial comparison of root-to-root loops {g) 
and paths getting closer to root by one level {h{) (see text). 

So we divide the expression for g{s) by the coefficient 
associated with the second step of g that hi does not 
take, 2a/-\/3, and multiply it by the coefficient of the 
step that hi takes instead, (1 -I- a)/-\/3. We also divide 
by the coefficient of the first step of g, not taken at all by 
hi, which is 1. In the special case s = 1 a single step is 
taken to the root from the first level, with a/-\/3. Finally, 
this path takes one step more as compared to g{s), so we 
use the expression for g{s + 1), starting from s = 3 since 



g{2) corresponds to the special case /ii(l). This gives us 
the expression for the ampUtude hi, 

l+g 
^^1(5) = -^Ms - 1) + -^^ 9{S + 1) X 1.6{3,5....}- 

Its transform is, using 1 + a + a^ — (as a = e^^*/'^), 
1 + a 



'^(^)=V!'+ 2a 



r ^^ — ^ 



t=3,5 
2 



(17) 
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The sum is the transform of g{t > 4) [304] , written as 
g — gt=2- With Eqs. (12), (16), and (17), the generating 
function for the amphtude of the process at the root is 



Hn{z) 



9iz)-V3z^ 



1 



(18) 



with g{z) given in Eq. (15). Now we need to invert this. 

B. Inverse transform: H„{t) asymptotic 

We take the inverse z-transform via an integral, and 
using Laurent expansion and the Residue theorem. 
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This integral is too complicated to yield a closed-form 
solution. Wc look for its asymptotic behavior in the form 



\z\=r 



1 — ^(z) Z* 



using the steepest descent method. The calculation is 
discussed in Appendix A. The asymptotic of the ampli- 
tude of the process at the root, starting from a level n in 
the tree, with r = t — n, is 






(20) 
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-i AiJ 
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Constants ci„ and C2n are linear in n, while 7 and A's are 
real constants (Appendix A). The probability is 

C2 - 2 Re (ciC^ g-i[27n+(Ai-A.)5] | 

l^"WI' ^ 4.^ 2» r3 ^' (2^) 



where C^ = |cip -|- |c2p ^ n^. The oscillations of the 
exponential term are rapid at large times (and/or n). 



and this function behaves as 



--3 



The n dependence 



is ~ n^/2" at large times. Also, we see that the walk 
is transient, in the sense introduced in [39, 40], since 
Er \Hnit)\^ is finite. 

The observed power law decay differs sharply from 
the exponential tail of the classical walk (Appendix B). 
On finite graphs built with binary trees, this exponen- 
tially slower decay may lead to algorithms with signifi- 
cant speedups. The treatment in this section is meant to 
lay the groundwork for such investigations. 

This behavior should also have general implications for 
physics of systems modeled with a binary tree. 



C. Inverse transform: Hn{t) computed 

The transform is defined as Hn{z) = J2t>n^n{t)z'^, 
and using its Taylor expansion and equating coefficients. 



Hnit) 



Hi'Hz)U=o 

t\ 



(22) 



This can be evaluated efficiently, for a range of values of t 
for a fixed n, providing the full amplitude. Symbolic cal- 
culation of derivatives (with Mathematica) allows for 
values of n in the thousands. We note that the amplitude 
shows an interference pattern, with the main peak fol- 
lowed by (much) smaller, rapidly diminishing, secondary 
peaks. Probability at the root with time is shown in 
Fig. 5, for n = 50. The shape does not depend on n. 
At long times this exact result can be compared with 
asymptotics (21), see inset in Fig. 5. The tail exhibits 
t^^ dependence, in agreement with Eq. (21). 



n=50 



numerical 



1.5 



asymptotic 




0.5 



time steps 

FIG. 5. (Color online) Probability at the root with time 
(n — 50), computed via Eq. (22). Inset compares this (points) 
with the steepest descent asymptotic (21) at same time steps 
(red squares). 



D. Comparison and algorithmic aspects 

We considered an algorithm of finding the root start- 
ing at the n-th generation in the tree. Here we compare 
the quantum and classical walks, and their estimated 
run times, using data calculated via Eqs. (22) and (B3). 
Probability peaks and times at which they are reached, 
for quantum and classical walks, are shown in Table I 
for a few initial levels n. Probability at the root for the 
classical walk is shown in Fig. 6. 




100 150 200 

time [steps] 

FIG. 6. (Color online) Classical walk, probability at the root 
(Appendix B). 

The best run time for a given n is estimated as fol- 
lows. The inverse of the probability is the average num- 
ber of times needed to run in order to hit the root; mul- 
tiplying it by its time gives the total running time to 
hit the root. We need its minimum over all time steps, 
min |i/|7J„(i)P}, usually the values for the peak. We 
use data up to n = 5000 in steps of 100 for the quantum 
walk, and up to n = 2000 in steps of 50 for the classical 
walk (numerical integration is more demanding). Using 
smaller increments in n does not affect results. We fit 
the natural logarithm (In) of run times with a polyno- 
mial and linear (Fig. 7). 

The polynomial fit establishes a linear trend (reached 
at n ~ a few hundred) , so the run time of the quantum 
walk is exponential in the initial distance from the root, 
~ e''". The same holds for the classical walk, and then 
the ratio of slopes of their linear fits compares their run 
times. 

This ratio does not change much over the whole range 
of data, being within a few percent of 2/3. Still, since an 
exponential complexity is fully felt at large n, the later 
portions of data are more relevant for algorithmic com- 
parison. For the last quarter of data ranges, indicated in 
Fig. 7 by lines fit through data, the ratio of quantum to 
classical slopes is « 0.685, within 3% of 2/3. Thus it ap- 
pears reasonable to conjecture the algorithmic speedup 
of the order of 2/3. (For a run time T for the classical 
walk, one expects the run time on the order of T^'^ for 
the quantum walk.) 

We note the behavior of peaks times with n. For the 
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polynomial fit, ~ n'' 
p, = 0.99907 ± 0.00003 (qw) ^ 
p, = 0.99875 ± 0.00002 (cl) 



X...-- 



linear fit, b x n 
bg = 0.475234 ± 0.000005 (qw) 
- b, = 0.69401 ± 0.00002 (cl) 
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FIG. 7. (Color online) Natural log (In) of run time witii 
initial level in the tree, for quantum (points) and classical 
(blue squares) walks. For plot clarity not all points are shown. 
Parameter values for polynomial (~ n^) and linear (bn) fits 
are shown, with lines through data used for fits (see text). The 
run time is ~ e**", and bg/bc — )■ 2/3 estimates the quantum 
over classical walk speedup. 



classical walk, t^ax(") = 
quantum walk i^x(") ~ 
fit of ^ 1.571, and for an 



3n — 8 (exact), while for the 
1.46n (where data allows for a 
-^ In n correction) . 



IV. 



CONCLUDING REMARKS 



Quantum walks are quantum processes with a specific 
mixing of states; particular unitary processes. In this 
vein, we propose to approach and study them using ideas 
from classical random walks with memory. 

We have demonstrated how a general framework for 
discrete-time quantum walks arises as a natural analog of 
a specific representation of a classical memory-2 Markov 
chain. Walks are implemented by constructing a local 
operator with no restrictions other than unitarity. The 
framework needs no "coin" degrees of freedom, is flexible, 
and is applicable to general graphs. This approach may 
make it easier to obtain walks on structures for which 
significant speedups are expected. 

The evolution operator works separately on each com- 
ponent of the amplitude, reducing the state space, and 
effectively deconstructing the amplitude [as in in the bi- 
nary tree example, Eqs. (7-9)]. This makes quantum- 
mechanical correlations and interference transparent in 
quantum walks, making their explicit study easier. It 
should also aid the use of quantum walks as a general 
tool for exploration and modeling of physical systems. 

In Sec. Ill we use the framework to build a symmet- 
ric discrete-time quantum walk on a semi-infinite binary 
tree. We start the walk at a level n in the tree and find its 
amplitude at the root as a function of time and n. The 
construction of the walk is simple, but the calculation 
is complicated by a (refiective) boundary. The gener- 
ating function of the amplitude is found explicitly, and 
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initial level n 


10 


20 


50 


100 


200 


500 


max j//„p (at t) 
max pt{n) (at t) 


6.8 X 10"* (16) 
1.2 X 10"'' (22) 


3.9 X 10"* (30) 
7.2 X 10"** (52) 


1.7 X 10"'^ (76) 
4.2 X 10"^'^ (142) 


5.3 X lO-'^* (150) 
2.6 X 10"^^ (292) 


7.6 X 10"^* (298) 
1.4 X 10"'^^ (592) 


5.1 X 10"'"* (738) 
4.5 X 10"'^^ (1492) 



TABLE I. Probability peaks at the root with their times, for quantum [|i/„(f)p] and classical [pt(n,0)] walks, with n. 



its asymptotic is found via the steepest descent method. 
The full solution is computed numerically. 

These solutions show interesting features. The asymp- 
totic decays in time by the power law (as opposed to the 
exponential tail of the classical walk), representing long- 
range correlations. This hints at significant speedups on 
restricted structures. The amplitude exhibits a damped 
interference pattern, with a distinct and sharp peak. 
In comparison with the classical walk, the probability 
peak is reached more quickly, and is orders of magnitude 
greater, already at small n. The run time for hitting the 
root on a semi-infinite binary tree is exponential with 
n for the quantum walk, as it is for the classical walk. 
It is still clearly slower in n and, following suggestive 
data trends, we conjecture the polynomial algorithmic 
speedup of the order of 2/3 over a classical walk. 



Appendix A: Steepest descent calculation summary 

Integrals suitable for analysis by the steepest descent 
method are typically of the form [41] 



I{k) = / f{uj) e'=*(") dw. 
Jc 



(Al) 



We use Eq. (19), where the exponent will be formed from 
powers of z. As z is always squared we first change vari- 
ables via z'^ = £,. Accounting for the double winding. 
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We now use g{£,)/^ = i^, to transfer some of the inte- 
grand's complexity into the exponent. 
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S, = (p{uj) = av3 - 
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Carrying out the substitution, we have 

{u> - \/3)" ipf 



I{t-n) 



if 2 (Jl^J^ 



(A2) 



in the form (Al), with tp ~ = e^"'"'^'^ \ and 

t — ft 
■hup, k= . (A3) 
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Keeping fi /(p will be useful. Consider the critical points 
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'73' 73" 



The logarithm's branch 



simple poles are at 

point is at ip 

what the contour must enclose, we can take any conve- 
nient branch. Two simple saddle points are 



1, and since this is not at LOp — l/\/3, 



(ln(p)/ = 



_ ^^ t V 2 _ ±i arctan \/2 



The main contribution to this integral comes from sad- 
dle points. A branch of the original integration contour 
\^/p\ = r can be chosen (via r) for use with steepest 
descent paths. There are no issues with deforming the 
contour, as no critical points are in the way, any branch 
of the logarithm is good, and k = ^^ G Z (Fig. 8). 
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FIG. 8. (Color online) Original integration contour (solid), 
steepest descent paths (dashed and dotted), and critical 
points. 

In the steepest descent method decreasing orders of 
contribution are computed mostly via expansion around 
saddles. (There are theorems and formulas for the first- 
order contribution, but here it is zero.) Around Ws we 
have <^{ijj) = ^{uJs) + If [*f(w)][^ (w - i:^^)^ -I- o{uj'^), and 
the usual change of variables $(a;) — ^(ws) = —y gives 



Tb^/y- 



(A4) 



where y is zero at the saddle and real along the steep- 



A pole of order *"^+^ is at (^ = ^ Wp = l/\/3. Two est descent path. We used ^ (ln(/3-i)" = i^J(\np) . 
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Now y = In [ip{uj)/ip{ujs)\ and dy = [(p{ujy /ip{uj)\ dui, 
and with (p{uj) = (p{uJs)e'^ the integral (A2) and (A3) 
becomes 



m 



,fe*(<^e) 



[uj{y)-V3]" 
1 -uj{y)ip{ujs)e.y 



'ky 



dy. (A5) 



Now we can directly expand around y = (to any order) , 
restricting integration to the line along the steepest de- 
scent path, close to saddles. The signs in Eq. (A4) cor- 
respond to the opposite directions from the saddle; we 
label +/— as "R/L." Substituting uj{y) and expanding, 



l^Aip; 



A 



W,s 



/ {l±B^)e'''ydy, 5^0(1), (A6) 
Jo 



1 - UJsfs 






For compactness we use (ps = 'fii^s) — e*^=, with Asi = 
arctan [(973+872)723], Xs2 — arctan [(973-872)723] — tt. 
Note that A{n) - (\/2)", as cj^ - v% = \/2e^'^% with 
7si7s2 — T [arctan(l/72) — tt] = =f(7 — tt), and we will 
extract tt later. The integral is dominated around lJs 
{6 m 0), so it can be formally extended, (5 — > 00, and we 
get / - A J^ (1 ± By/y) e-'^y dy. This results in 
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Subtracting contributions along opposite directions, the 
first non-zero order for either saddle is 



/, 



(V2)", 



,* "fsfl' 



{as + ds n) 



-i Xgk 



kVk 
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We broke up the A^Bn term found in Ir~ II, to show 
the structure of n dependence, where a., — ,-, '^f A2 and 



b7^ 



"'* (l-WsVs)('^«-V3) 

TT from 7s, and will use e^^- 



are of order 



1. Here we extract 



= (—1)6^*'''. Contributions 



for saddles are subtracted (for consistency of ±y^ direc- 
tions) and, with Cs^„ = a^ -I- dgn, we get Eq. (20). 

The full expansion of the integral ( A5) results in nested 
sums of a Gamma function. This cannot capture the 
peak of the amplitude though, and is not needed for our 
asymptotic analysis, so we do not pursue it here. 



Appendix B: Random walk on a semi-inflnite line 

For completeness here we provide the application of the 
method developed in [42, 43] (based on Karlin-McGregor 
spectral approach to random walks) to a classical walk 
on a binary tree. 

If P is a reversible Markov chain over a sample space 
r2, and tt is a reversibility function (not necessarily a 
probability distribution), then P is a self-adjoint operator 
in ^^(tt), the space generated by the inner product 

< La >iT= X! fi^)9{x)TT{x) 
xes 



induced by tt. If P is tridiagonal operator (i.e. a nearest- 
neighbor random walk) on f2 = {0, 1,2,...}, then it must 
have a simple spectrum, and is diagonalizable via orthog- 
onal polynomials, as it was studied in the 1950s and 1960s 
by Karlin and McGregor. There the extended eigenfunc- 
tions Qj{\) {Qo = 1) are orthogonal polynomials with 
respect to a probability measure ip and 



Pti_i,j) 



A*Q,(A)Qj(A)rf^(A) Vi,jer!, 



where ttj (ttq = 1) is the reversibility measure of P. Con- 
sider the following Markov chain 
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Orthogonal polynomials are obtained via solving a simple 
linear recursion: Qo = 1; Qi = A, and 



)„(A)-ci(A)pl'(A) + c2(A)p^(A), 



(A7) where pi(A) = 



_ X + ^X^-ipq ^^^ _^^^ _ X-^X^-4pq 



2p 



and p2(A) = 



2p 



are 



the roots of the characteristic equation for the recursion. 



and ci 



and C2 



X-pi 



Now ttq = 1 and 



— p_ 



- {n > 1). Also, we observe that 

|P2(A)| > V^Tp on [-1, -2y^), 
\P2iX)\ < y/q/p on(2v^,l], 
|P2(A)| = y/q/p on [-2^^, 2^], 



and P1P2 — ^. The above will help us to identify the 
point mass locations in the measure ip since each point 
mass in ip occurs when J2k''^kQkW < °o- Thus we 
need to find all A e (2y/pg, 1] such that ci(A) — and 
all A G [— 1,— 2y/pg) such that C2(A) = 0. But there 
are no such roots, as ci(— 1) = and 02(1) = 0, while 
— 1 ^ {2^/pq,l] and 1 ^ [—l,—2y/pq). Thus there are 
no point mass atoms in ^, and the mass of ip must be 
continuously distributed inside [— 2y^, 2y^]. In order 
to find the density of V' inside [—2y/pq, 2^/pq] we need to 
find [eg, (P — s/)~^eo] for Im(s) 7^ 0, i.e. the upper left 
element in the resolvent of P. 

Let (ao(s),ai(s), . . . )"^ = (P — s/)^^eo, then 

— sao + ai = 1, and ga„_i — sUn + pun+i = 

Thus a„(s) = aiPi(s)" + a2P2(s)", with a, = 2[T)-Ms] 
anda2 = p7(°)-^pil^) • Si'^ce (ao,ai,...) e^2(C,7r), 



as n -> 




Hence when |pi(s)| 7^ |/02(s)|, either ai — ot a2 — 0, 
and therefore 



ao(s) 



''|pi(«)l<^ 
Pi(s) - 



P2(s) - 



(Bl) 



Also d'il}{z) — ip{z)dz, where (p{z) is an atom-less density 
function over [—2^^/pq,2y/pq], and 



ao(s) = 



+2\/M 



-2Jpq 



dTp{z) 



-2VM 



-2Jpq 



ip{z)dz 



Next we use the following basic property of Cauchy trans- 
forms Cf{s) = 2^ /h z-s^ *^^* *^^'-^ ^^ derived using the 
Cauchy integral formula, or similarly, an approximation 
to the identity formula: 



C4 



C_ 



(B2) 



Observe that the curve in the integral need not be in 
for C+ - C_ = / to hold. Here 

C+f{z) = lim Cf{s), and 

s^z: lin(s)>0 

C_/(z) = lim Cf{s), 

s^z: lin(s)<0 

for all z G M. The relation (B2) implies 



^p(x) 



27ri \ s=a;+ie: 
* e-i-0 



lim ao(s) — lim ao(s) 

-x^ie: s—x~ie: 
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for all X £ {—2^/pq^2^/pq). Recalling (Bl), we express 
which in turn simplifies to 



■\/4pg- 



(p(a;) = <^ 277g(l-x2) 

b 



if x e (-2^/pg,2^), 
otherwise. 



Here (^(x) always integrates to 1 over {—2y/pq^2^fpq). 
Now 



Pt{n,0) = 



-2vM 

-I-2VP9 



A*g„(A)^(A)dA 



A*(cip?-fC2p^) 



(pi - P2)d\ 
2m{pi - A)(p2 - A)' 



and therefore, since ci = -^^ — and C2 = 



A-pi 



Pt(",0) 



27ri 



+2VM 



2v^ 



P'i 



Pi 



/02 - A pi - A 



dA. (B3) 



This can be treated as a complex integral, for example, 
with steepest descent. But one can observe directly in 
Eq. (B3) that the tail of pt(n, 0) decays as {2y/pqy when 
t -^ +00. Thus, using p = 2/3 and q = 1/3 for the 
classical symmetric walk on the semi-infinite binary tree, 
the decay rate will be (2-\/2/3)*, giving us the exponen- 
tial asymptotics. The probability integral (B3) can be 
efficiently evaluated numerically (see Fig. 6). 
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